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ABSTRACT 

We investigate the effects of discontinuous mass loss in recurrent outburst events on the 
long-term evolution of cataclysmic variables (CVs). Similarly we consider the effects 
of frictional angular momentum loss (FAML), i.e. interaction of the expanding nova 
envelope with the secondary. The Bondi-Hoyle accretion model is used to parameterize 
FAML in terms of the expansion velocity Vexp of the nova envelope at the location of 
the secondary; we find that small Vexp causes strong FAML. 

Numerical calculations of CV evolution over a wide range of parameters demon- 
strate the equivalence of a discontinuous sequence of nova cycles and the corresponding 
mean evolution (replacing envelope ejection by a continuous wind), even close to mass 
transfer instability. A formal stability analysis of discontinuous mass transfer confirms 
this, independent of details of the FAML model. 

FAML is a consequential angular momentum loss which amplifies the mass transfer 
rate driven by systemic angular momentum losses such as magnetic braking. We show 
that for a given Wexp and white dwarf mass the amplification increases with secondary 
mass and is significant only close to the largest secondary mass consistent with mass 
transfer stability. The amplification factor is independent of the envelope mass ejected 
during the outburst, whereas the mass transfer amplitude induced by individual nova 
outbursts is proportional to it. 

In sequences calculated with nova model parameters taken from Prialnik & Kovetz 



( 1995 ) FAML amplification is negligible, but the outburst amplitude in systems below 
the period gap with a white dwarf mass ~ 0.6Mq is larger than a factor of 10. The mass 
transfer rate in such systems is smaller than 10~^^MQ/yr for ~ 0.5 Myr (~ 10% of the 
nova cycle) after the outburst. This offers an explanation for intrinsically unusually 
faint CVs below the period gap. 
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1 INTRODUCTION 

Cataclysmic variables (CVs) are short-period binary sys- 
tems in which a Roche-lobe filling low-mass main-sequence 
secondary transfers mass to a white dwarf (WD) primary. 
The transferred matter accretes onto the WD either through 
a disc or a stream and slowly builds up a hydrogen-rich sur- 
face layer on the WD. With continuing accretion the pres- 
sure at the bottom of this layer increases, and hydrogen 
burning eventually starts. The thermodynamic conditions 
at ignition determine how the burning proceeds (e.g. Fuji- 
moto 1982). If the degeneracy is very high, a thermonuclear 
runaway occurs, leading to a violent outburst terminated 



by the ejection of all or most of the accumulated envelope. 
Classical novae are thought to be objects undergoing such 
an outburst (cf. Livio 1994 for a recent review). Ignition at 
moderate or weak degeneracy causes strong or weak H shell 
flashes, whereas stable or stationary hydrogen burning re- 
quires fairly high accretion rates ^ W'^ Mq/jv which are 
not expected to occur in CVs. 

Mass transfer in CVs is driven by orbital angular mo- 
mentum losses which generally shrink the binary and main- 
tain the semi-detached state. The observed properties of 
short-period CVs below the CV period gap (orbital period 
P ^ 2 h) are consistent with gravitational wave radiation as 
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the only driving mechanism. A much stronger angular mo- 
mentum loss, usually assumed to be magnetic stellar wind 
braking, is needed for systems above the gap (P 3 h). The 
assumption that magnetic braking ceases to be effective once 
the secondary becomes fully convective in turn provides a 
natural explanation for the period gap as a period regime 
where the systems are detached and therefore unobservable 
(Spruit & Ritter 1983, Rappaport et al. 1983). The resulting 
typical mass transfer rate X in CVs is X ~ 5 x 10~^^MQ/yr 
below the gap and X ~ 10~^ — lO~*Af0/yr above the gap 
(see e.g. King 1988, Kolb 1996, for reviews). 

For negligible wind losses from the system the accretion 
rate is essentially the same as the transfer rate. With the 
above typical values the H ignition on the WD turns out to 
be degenerate enough to cause more or less violent outbursts 
(e.g. Prialnik & Kovetz 1995). As the mass to be accumu- 
lated before ignition is very small (AMig ~ 10"'^ - 10~^ Mq) 
the outbursts recur on a time tree = AMig/X much shorter 
than the mass transfer timescale which determines the long- 
term evolution. Studies of the secular evolution of CVs make 
use of this fact and replace a sequence of nova outbursts with 
given recurrence time tree and ejected envelope mass AA/ej 
by a continuous isotropic wind loss from the WD at a con- 
stant rate AMej/Uec- 

Such a procedure obviously neglects any effect that nova 
outbursts may have on the long-term evolution. These are 
in particular 

(i) The evolution of the system is not continuous but charac- 
terized by sudden changes of the orbital parameters, causing 
the mass transfer rate to fluctuate around the continuous 
wind average value (see Sect. 2). It is not a priori clear if 
the continuous wind average properly describes the system's 
evolution close to mass transfer instability. 

(ii) At visual maximum (and the following decline) nova 
envelopes have pseudo-photospheric radii of typical giants, 
i.e. much larger than the orbital separation. Therefore the 
secondary is engulfed in this envelope and possibly inter- 
acting with it. Drag forces on the secondary moving within 
the envelope can lead to frictional angular momentum loss 
(FAML) from the orbit, and accretion of envelope material 
onto the secondary could increase its photospheric metal 
abundances (pollution) . 

(iii) The H burning hot WD is extremely luminous (~ 
1O*L0, compared with accretion luminosities ~ li©) for 
as long as a few years. This might drive additional mass loss 
from the secondary star. 

In this paper we will focus on the effects of mass loss 
discontinuities and FAML. We neglect irradiation as it does 
not last long enough to influence the long-term evolution. 
Pollution is expected to be important only for metal poor 
secondaries (Stehle 1993); we neglect it altogether. 

In Sect. 2 we formally derive the continuous wind aver- 
age and consider the stability of mass transfer in the pres- 
ence of nova discontinuities analytically, with FAML of ar- 
bitrary strength. We review previous studies on FAML and 
follow Livio et al. (1991) to derive a simple quantitative 
model for FAML in Sect. 3. Using this description we per- 
form numerical calculations of the long-term evolution of 
CVs with various strengths of FAML, both for sequences 
of nova outbursts and the continuous wind average. Re- 
sults of such computations where the FAML strength and 
the ejected mass per outburst have been varied systemati- 



cally are shown in Sect. 4. Sequences with FAML parameters 
taken fr om th e consistent set of nova models by Prialnik & 
Kovetz (1995) are shown at the end of Sect. 4. Section 5 
discusses our results. 



2 MASS TRANSFER STABILITY AND 
CLASSICAL NOVAE 

We begin by investigating how mass transfer discontinuities 
induced by nova outbursts affect mass transfer stability. In- 
troducting the well-known conservative mass transfer sta- 
bility criterion we develop a formalism to extend its appli- 
cability to the discontinuous case. The strength of FAML 
enters as a free parameter. 



2.1 Conservative and continuous CV evolution 



Following Ritter (1988) the mass transfer rate X in a CV 
can be approximated by 



X = -M2 = Mo exp 



R2-RR 



(1) 



(note that X is always positive). Here both Mq ~ 
10"** M© yr"^ and the ratio e = Hp/R2 ~ Hp/Rn ~ 10"* of 
the photospheric pressure scale height Hp to the donor's ra- 
dius J?2 are roughly constant for the secondaries under con- 
sideration. The secondary's Roche radius 7?r can be written 
as a fraction /2 of the orbital distance a, 



Rr = /2(g) a- 



(2) 



/2 depends only on the mass ratio q = M1/M2 and is for 
g > 1.25 given by /2(g) ~ (8/81(1 + q)Y^^ to better than 
2% (Paczyhski 1971). 

To assess how X changes with time and to find a sta- 
tionary value for X (where X = Q) vie have to consider both 
R2 and -Rr. From the total orbital angular momentum J, 



and the derivative /32 
find with (0) 



din /a /ding (= 



(3) 

-g/3(l + g)) we 



(4) 



where Mi, M2 and M = Mi -I-M2 denote the primary mass, 
donor mass and total mass, respectively. In (Uh we formally 
separate "mass-loss related" angular momentum loss Jm, 
caused by mass which leaves the binary and carries a certain 
specific angular momentum j = i'{J /M) (quantified by the 
dimensionless parameter v), i.e. 



Jn: 



-jM 



M 



(5) 



and "systemic" angular momentum loss Jsys which oper- 
ates without (noticeable) mass loss, e.g. gravitational wave 
radiation and magnetic braking. By combining the mass- 
changing terms (U|) is usually rewritten as 



Rr . M2 , „ J.ys 

r^-^^'m;^^—' 



(6) 
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Figure 1. Reaction of the secondary's Roche radius to mass 
transfer and mass loss, as a function of M2 for fixed Mi = 1.2Mq . 
Plotted is the mass-radius exponent for conservative mass trans- 
fer {C^ - full line), isotropic wind with r] = 0, u = l/q (fp^ from 
(^; long dashed), and isotropic wind with rj = 0, u = q (short 
dashed). Also shown is for various values of (see Sect. 3.1). 



thus defining tlie mass-radius exponent ("r of tlie sec- 
ondary's Roclie radius. In the simple case of conservative 
mass transfer where M = constant, i.e. M2 = —Mi, we find 
from (|i) 



(7) 



In the more general case of an isotropic wind loss at a rate 
M = (1 - J7)M2, i.e. Ml = -r)M2, we obtain 



Cr. = (1-.)^+.--2-/?Ji + 2 
l+q q \ Q 



(8) 



an expression especially useful when ri and v are constant. 
Figure ^ depicts and ("r with 77 = and either v — l/q 
ox V = q as a. function of M2 for Mi — 1.2Mq. 

Similarly, it is standard practice to decompose 
the secondary's radius change into the adiabatic re- 
sponse (adM2/M2 and the thermal relaxation K = 
{d\nR2/dt) 

Af — const. ; 



(9) 



Fig. y shows the adiabatic mass-radius exponent (^ad as a 
function of stellar mass for low-mass ZAMS secondaries 
(Hjellming 1989). Fully convective stars and stars with deep 
convective envelopes (M2 ^ O.SM©) have C,aA — —1/3. 

Using these definitions the time derivative of X accord- 
ing to (|l|) becomes 



X = X [A-BX] 
with 



(10) 



Figure 2. Reaction of the secondary's radius to mass loss as a 
function of stellar mass M2 - Full line: adiabatic mass-radius index 
Cadi dotted line: thermal mass-radius index C,e- 



and 



^? = -ri- (Cad - Cr) 



(11) 



(12) 



The binary attempts to settle at the stationary mass 
transfer rate 



Xs^ jr =M2 — 

B Cad — Cr 



(13) 



{X = for X — Xs). The stationary rate is stable if 
dX/dX < at X = Xs, i.e. if the system opposes any in- 
stantaneous perturbation in X. This translates into Z3 > 0, 
hence the familiar criterion 

Cad — Cr > <4> mass transfer (dynamically) stable . (14) 

A similar stability criterion, Cc — Cr > Oi against thermal 
timescale mass transfer can be derived (see e.g. Ritter 1996), 
where Ce is the thermal equilibrium mass-radius exponent 
(~ 0.85 for low-mass main-sequence stars, see Fig. 



2.2 Non— stationary mass transfer 

The differential equation (^) governs the evolution of the 
mass transfer rate as a function of time t. For constant A 
and B the solution of ( p^ ) is 

where Xo = X(t = 0) denotes the initial value of X. Hence 
any changes of the transfer rate proceed on the character- 
istic timescale tc = l/A ~ 0{etj), i.e. is a small frac- 
tion e of the systemic angular momentum loss timescale 
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In reality both A and B change with time, typically on 
a the secular timescale tj, but for most practical cases where 
we investigate mass transfer stability against short timescale 
perturbations A and B can be considered as constant. 

Formally (1^) is a solution of (|l^) even if ^ < (i.e. 
\K\ > 2 Jays/ J ). In this case, assuming dynamical stability 
(S > O), Xs is negative and no longer a stationary value for 
the mass transfer rate. Rather (|l^) shows that in this case 
X decreases exponentially, i.e. the system detaches (note 
that Xo is a physical mass transfer rate and as such always 
positive). If finally the mass transfer is unstable {B < 0) then 
we see from ( |lo| ) that unless A has a large negative value 
the transfer rate grows. The growth timescale is initially tc, 
but becomes shorter and shorter with further increasing X. 



2.3 Nova— induced discontinuities 

So far we considered only continuous (though not neces- 
sarily stationary) mass transfer. Events that change orbital 
parameters on a timescale much shorter than the character- 
istic time etj ~ 10'' yr for reestablishing the local stationary 
value Xs may be regarded as discontinuous and instanta- 
neous. Nova outbursts certainly belong to this category; the 
nova envelope expands beyond the orbit within less than a 
few 10^ days after ignition and returns within 1 — 10 years 
(actually these are upper limits for the slowest novae). In 
the following we apply the above formalism separately to 
the outburst phase and the inter-outburst phase, then com- 
bine them to describe the full nova cycle. 



2.3.1 Outburst phase 

As a result of the outburst the envelope mass AAfcj is 
ejected. We expect that the specific angular momentum car- 
ried away by AMoj is higher than the specific orbital angu- 
lar momentum of the WD due to dynamical friction of the 
secondary orbiting within the envelope (FAML). Hence we 
write for the discontinuous change of the orbital angular 
momentum 

AJ=-(l/l+l/FAML)-^AMoj, (16) 

where = M2/M1 accounts for the specific orbital angular 
momentum j = ji of the WD (represented as a point mass). 
The free parameter ^taml measures the strength of FAML 
and will be estimated in terms of a simple model for the 
frictional processes in Sect. 3 below. 

Correspondingly, using (|) with AM = AMi = -AAfoj 
and AM2 = (i.e. neglecting the small amount of mass 
accreted onto the secondary as well as any mass transfer 
during this short phase), the change A_Rr — i?R(post) — 
_RR(pre) of the Roche radius is 



( ARr 



AMei 



Ml \l + q 



-P2 



_2q_ 

1 + g 



i^FAML 



(17) 



2.3.2 Inter- outburst phase 

After the outburst mass transfer continues. Once the mass 
accreted onto the WD exceeds the critical ignition mass 
AMig > the next outburst occurs. For simplicity we as- 
sume conservative mass transfer during the inter-outburst 



phase. According to (|6|) the total change of the Roche radius 
in the inter-outburst phase is 

/ A^\ ^ ^AA/ig 

V / inter M2 



n A(/sys 



(18) 



where AJgyg is the total change of the orbital a.m. due to 
systemic losses between outbursts. 

2.3.3 Combined description 

For the study of the long-term evolution of CVs it is conve- 
nient to replace the sequence of nova cycles with mass trans- 
fer discontinuities by a mean evolution where the mass AAfej 
is regarded as being lost continuously at a constant rate 
over the cycle in form of an isotropic stellar wind carrying 
the specific orbital angular momentum (vi + vpAMh) J/M. 
If Cr denotes the corresponding Roche lobe index then we 
have from (^) 

/Ai?R- 



_ p AAfa 

V i?R y 

total 



r, AJsys 

2 — 



J 



(19) 



again with AM2 = — AAfig because we neglect any change 
of A/2 during outburst. As the total change of the Roche 
radius over a complete cycle is A_Rr (total) — A7?R(out) -|- 
A_Rr (inter), comparison with ( p^ ) and ( p^ gives 

- + -AM,JM2 ■ 
This can be written as 



Cr 



q 

l-r," 



2q 

1 + g 



I (21) 



— 2 UpAMh 



Here 1 - ry" = AM1/AM2 = AA'/oj/AA/ig specifies the 
change of the WD mass during outburst in units of the mass 
lost from M2 during the inter-outburst phase. Equation ( |2l| ) 
in fact is equivalent to (^) with rj = rj" and = i^i + i'faml ; 
77" determines if the WD mass will grow or shrink in the 
long-term evolution. 

The generalization of ( ^ ) which allows for non- 
conservative mass transfer between outbursts and a change 
of AI2 during outburst is given in the Appendix. 

2.4 Stability of discontinuous mass transfer 

As was shown in Sect. 2.1 mass transfer in the fictions mean 
evolution which mimics the effect of nova outbursts by a 
continuous wind loss is dynamically stable if 



Cad 



Cr > . 



(22) 



Clearly, in the case of discontinuous mass transfer with nova 
cycles the stability considerations leading to this criterion 
are no longer applicable, as the system never settles at a 
stationary transfer rate. After an outburst the transfer rate 
X follows the solution (^5|) and approaches the conservative 
stationary rate X^ = A/B'^ (where B'^ is B according to ( |l2| ) 
with Cr)- At the next outburst X changes discontinuously 
from the pre-outburst value Xprc to the (new) post-outburst 
value Xpost, i.e. increases (or drops) by a factor 



Xp 



: exp 



1 fAR 



( Rr 



(23) 
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see Eqs. (^) and (|l^)- In the following we consider mass 
transfer stability in the presence of nova cycles and derive 
a generalized stability criterion that reduces indeed to the 
simple form (|2^). 

In our approach we compare the mass transfer rate X^^^ 
and Xpf^ immediately before subsequent outbursts i and i+1 
in a sequence of cycles with constant AMig, AMej, and 
A. Eq. ( p^ with Xq — X^^^^ describes X as a function of 
time between outbursts. At time t = tree when the outburst 
i + 1 ignites, the mass transferred since the last outburst is 
just AMig, i.e. 



X{t) dt = AMig 



(24) 



This can be solved for the outburst recurrence time t^ec (cy- 
cle time) 



1 

A 



In 



cxp {B'= AMig 



(25) 



Inserting tree for t in (|l3) gives the new pre-outburst value 

xt+^ ■ ■ 



X, 



pre — ^(^rcc) , 

Ito = X^ [1 - exp(-B=AMig)l 
+X;,eexp(-B^AMig) , 



(26) 



where we used ( |23| ) to replace X^^^^ by X^^cj ( P0| ) and the 
definition of analogous to B'^, i.e. B^ — (^ad — CK)/^^2- 
Hence the relative change of the pre-outburst value is 



:f{x') = 
x_^ 

X 



X 



i+1 



X' 



(27) 



I [1 - exp(-S''AAfig)] - [l - exp (-B^AMg 

A_ 

X^ 



AMig 



B' 



(the subscripts for all X have been dropped as they all read 
"pre"). We will refer to as the growth function. 

For > the mass transfer rate grows from outburst 
to outburst, for < it decreases, on a timescale t^acIT . 
Formally, .^^ = for X = X, with 



X = 



1 - exp (-S'^AAfig) 
1 - exp (-ZS^'AA/ig) 



A_ 

BP 



(28) 



X represents the stationary pre-outburst mass transfer rate 
if it is positive, i.e. if A and B^ have the same sign. For 
it to be stable we require dTjdX^ < at X' = X so 
that the system evolves back to X after a perturbation in- 
creases/decreases X. This is the case if, and only if, ^ > 0. 
Then the stationary solution exists only if B^ > 0, in all 
other cases the system either detaches, or the transfer rate 
grows unlimited. Figure |^ summarizes schematically the for- 
mal functional dependency of on X' for various com- 
binations of the signs of A and (see Table 0). Hence 
the generalized stability criterion for discontinuous nova- 
amplified mass transfer is B^ > 0, which is indeed equiv- 
alent to (I22I). Remarkably, this condition is independent of 
the sign of B'^, i.e. the system can evolve in a stable man- 
ner with formally dynamically unstable conservative mass 
transfer between subsequent outbursts. 

The analysis of the system's behaviour once it violates 
condition (|22[) is not straightforward. The main problem is 




log ( XV Mo yr-' ) 



Figure 3. Growth function JP, defined in Eq. (|27[), as a func- 
tion of pre— outburst mass transfer rate X* (schematically). See 
Table |^ for a description of the different curves. Only the solid 
curve represents a stable stationary solution X (marked by a filled 
diamond), all other cases lead to a runaway or a turn-off. 



Table 1. Branches of possible solutions of Eq. (^) as shown in 
Fig. I 



LINE 


SIGN OF 


STATIONARY 


STYLE 


A 


SOLUTION 


full 


+ + 


yes (stable) 


long dashes 




yes (unstable) 


dotted 


+ 


no 


short dashes 


+ 


no 



that A is expected to change its sign (from positive to neg- 
ative) at about the same time as B^ . This can be seen from 
an estimate of the thermal relaxation term K if the real 
evoluti on is replaced by the continuous wind average: Stehle 
et al. ( 1996| ) have shown that the secular evolution rapidly 
converges to an attracting evolutionary path characterized 
by 



R2 _ M2 
R2 ^ A/2 ' 

whatever the initial co nfig uration of the s; 
stationarity, equating ( p9| ) with (^, and 
after elimination of M2/M2 



A: 



2 Jsys Cai 



SR 



(29) 

stem. Assuming 
with (ph gives 



(30) 



as an estimate for A. Hence ^ > during a phase of stable, 
stationary mass transfer, but ^ < once ("ad < Cr- 
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In practice this means that when a stable system, char- 
acterized by the full line in Fig. ^, approaches instability, the 
curve becomes flatter and flatter with X almost constant, 
and then inverts to the (unstable) long-dashed line. 



3 A SIMPLE MODEL FOR FAML 

In the previous section we derived analytic expressions de- 
scribing stationary mass transfer and mass transfer stabil- 
ity in the presence of nova cycles. A non-zero FAML effect 
was allowed, and the strength of FAML was treated as a 
free parameter j^aml. To proceed further and underline the 
above findings with numerical examples we quantify vpAMh 
by adopting a model for the frictional processes leading to 
FAML. 



MacDonald (198C) discussed for the first time the possi- 



bility that an interaction between the secondary and the ex- 
tended envelope might tap energy from the orbit and thereby 
reduce the nova decay time. He determined the angular mo- 
mentum transferred from the orbit to the envelope by de- 
scribing the nova envelope as a polytrope with index 3 at 
rest, and restricting the combined nuclear, internal and fric- 
tional energy generation to t he Ed dington luminosity (Mac- 
Donald 1986). Shara et al. ( |l986| ) and Livio et al. ( il99l| ) 
used a more direct way, based on Bondi-Hoyle accretion, to 
estimate the transfer of angular momentum. Their approach 
(which we adopt below) offers a very simple analytic treat- 
ment, as the explicit expression for FAML has only one free 
parameter, the envelope expansion velocity iicxp, measur- 
i ng th e strength of FAML. More recently, Kato & Hachisu 
(1994) have again confirmed their earlier result that FAML 
has only minor infiuence on the decay phase of nova out- 
bursts. This can be understood in terms of the high ex- 
pansion velocities in their models and will be discussed in 
Sect. 5. On the other hand, Lloyd et al. (1997) showed that 
common enevlope evolution can contribute to the shaping 
of the nova remnant. 



3.1 



FA ML de scription according to Livio et 
al. (|l99l|) 



In Bondi-Hoyle accretion linear momentum is transferred 
through the drag force 



-Fdrag = — Cdrag(A^, 7gas) - TT R\ /5 UrolWrcl , 



(31) 



medium. The dimensionless drag coefficient Cdrag 
Mach number M _and specific heat ratio 7gas 



where v^ci is the gas stream velocity at infinity, _Ra the ac- 
cretion radius (see below) and p the density of the accreted 

varies with 
and is of order 

unity. Equation (^) is strictly valid for the highly super- 
sonic case, but also applicable for TVf ^ 1 if an adequate 
interpolation for the accretion radius is used, e.g. 

2GM2 



Ra = 



(Shima et al. 1985); cs is the local sound speed. Detailed 
3-dirn . hydrodynamical calculations summarized in Ruffert 
( 1995 ) confirm that (|3^) describes the subsonic case as well 
when nR\ is replaced by the geometrical cross-section. But 
these simulations also reveal additional complexities not rep- 
resented in the simple form (pll), e.g. the influence of the 



(32) 




exp 



rel 



Figure 4. Sketch of the simplified geometry assumed to derive 
i^FAML- The WD (Ml) is at rest and orbited by the secondary 
(M2) at distance a with velocity fscc- The envelope expands 
spherically symmetric from the WD, and has constant expansion 
velocity Wexp at the location of the secondary. 



accretor's size, even in the supersonic case. Most important, 
the dr ag coe fficient can deviate significantly from unity. Kley 
et al. (1995) point out that radiation pressure may reduce 
Cdrag by a factor of ~ 20 because it increases the sound speed 



and thus lowers M to the subsonic case. 

Given these uncertainties, and the fact that idealized 
Bondi-Hoyle accretion is merely a rough approximation to 
the situation of a secondary orbiting in an expanding nova 
envelope, the quantitative expressions for FAML derived be- 
low are order of magnitude estimates only. 

We now consider the FAML situation in a simplified 
geometry (Fig. the WD is at rest in the centre of a 
spherically symmetrical expanding nova envelope. Then the 
secondary's orbital velocity 



G (Ml + M2) 



(33) 



(typically Vscc — 400 — 500 km/s) is perpendicular to the 
wind expanding with ZJexp- Hence the velocity of the accre- 
tion flow with respect to the secondary is 



cl — V ^scc 



(34) 



As M2/MQ ~ R2/RQ the ratio of accretion radius (^ to 
stellar radius R2 is 



Ra ^ f 360km/s 



(35) 



(where v is the quadratic mean of iigcc, Vcxp and cs), sug- 
gesting that Ra < R2 for most cases. 

Setting Cdrag = 2 and Ra ~ R2 we obtain for the loss 
of angular momentum due to friction from the tangential 
component of the drag force ( ^l| ) 

jpAML = —\a X Fdragl = —a n R^ p Vscc Urol- (36) 

Using the continuity equation 

- Ml = 47r r V Wexp (37) 
to eliminate p, (^) and (^) then give 
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Jfaml = - [/2(g)]' — Ml VGM^. 

4 Vcxu 



(38) 



If most of the mass and angular momentum can be as- 
sumed to be lost at roughly constant rate and velocity (Livio 
at al. 1991) we finally obtain the total FAML, integrated over 
the outburst, 



A Jfaml 
where 



- -J ^ A{q) B{K,), 



B{K,) = 



Vrcl 



and 

Vscc 



J^FAML 



Comparison with ( |l^) then shows that 

1 + g 



'-Aiq)B{K, 



i.e. the mean specific a.m. jej of the ejected material is 

Jcj=(i+i^FAML)^. 



(39) 

(40) 
(41) 

(42) 
(43) 
(44) 



The functional dependence on Ki in Eq. ( pl[ ) , shown in 
Fig. ^(top panel), leads to the unphysical limit — A Jfaml — + 
00 for a very slow envelope expansion Uoxp 0. This is 
mainly due to the neglect of the envelope spin-up which 
would effectively reduce the relative velocity (Livio et al. 
1991). 

The secondary cannot spin the envelope up to speeds 
faster than corotation, so jcj is certainly smaller than ji + 
a^u! (where uj is the binary's angular velocity). If only a 
certain fraction of the envelope is spun-up to corotation 
this upper limit is correspondingly smaller, e.g. 

jmax = ji + {R2/a)a^Lj (45) 

if only the torus traced by the secondary's orbital motion 
corotates. 

Equation (^) specifies our model of FAML in terms of 
a single parameter, KiJot Vcxp) and determines the value 
of according to Eq. (£^). We use (^5|) as a physical upper 
limit on the effect of FAML. 
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Figure 5. Top panel: the function B(Ki) defined in (^). Lower 
panels: the specific angular momentum of the ejected matter as 
a function of i (or "Uoxp , obtained by assuming a constant Vbcc = 
400 km/sec), measuring the strength of FAML, for different mass 
ratios q = M1/M2, according to (k4). The dashed lines mark the 
specific a.m. ji and j2 of the WD and the secondary, respectively. 



3.2 The effect of individual nova outbursts 

To illustrate the effect of FAML in the parametrization ( |43| ) 
we consider the resulting relative change Ax/x — (a;(post) — 
x(pre))/a; of a system parameter 2: as a result of the nova 
outburst. For the orbital distance we obtain, similar to (171), 



Ml \l + q 
10"* (1-0.2 5(^-1)) 



AMoj/10-*Mc 



nil 



which translates into 



AA/cj I 2_S_ ^ s A{q) B{Ki 



V P /out Ml \ 1 + q 
10~* (2 - 0.6 B(J^i)) 



AMoj/10-*Me 



nil 
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(46) 



(47) 



for the orbital period P. The approximate expressions in 
(H) and (|^) make use of the weak dependence on q; mi is 
Mi/Mq. For a given ejection mass and white dwarf mass, 
B{Ki) critically determines the magnitude of the orbital 
change. From (^^ and ( |l7[ ) the change of the mass transfer 
rate X = -M2 is 



f — ) 



(48) 



AMcj/lO-*Af0 



(1 - 0.2 B{Ki)) ■ 



mi 



where we used e ~ 10 * to obtain the last line. Figure 
illustrates Eq. (H). 
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Fi gur e 6. Outburst amplitude of the mass transfer rate according 
to (ha) for various strengths of FAML (as labelled) and AM^j = 
10-^Mq, Ml = 1.2 Mq, Hp/-Rr = 1.4 10-4. 



Depending on the strength of FAML and the mass ra- 
tio the changes in system parameters can be either positive 
or negative. The critical -ft^i(g) where the amplitudes vanish 
are shown in Fig. ^. The critical lines for P and X almost 
coincide and separate the two different outburst types, those 
which increase the mass transfer rate and decrease the or- 
bital period, and those which decrease the transfer rate and 
increase the period. 



3.3 The effect on the long— term continuous wind 
average 

To consider the effect FAML has on the continuous wind 
average evolution we first note that FAML is a consequential 
angular momentum loss (CAML; cf. King & Kolb 1995) as 
jpAML is proportional to the mass transfer rate, see ( |39[ ) 
with AMcj = XpAMLtroc (whcre X denotes the continuous 
wind average mass transfer rate). As such FAML amplifies 
the transfer rate which would be driven by systemic angular 
momentum losses alone in the absence of FAML. 

An analytic estimate of this amplification factor, the 
ratio Xfaml/A' of the mass transfer rate with and without 
FAML, can be obtained if we assume that the system fol- 
lows the corre spond ing uniform evolutionary track found by 
Stehle et al. ( [l996| ) also in the case with FAML. This is a 
reasonable assumption as long as the amplification is mod- 
erate and the resulting mass transfer timescale longer than 
the thermal time of the secondary. The radius reaction of 
the secondary along its track is given by (^^. Hence using 
Rr/Rr = R2/R2 and (El) we have 



^FAML (Cc ~ Cr) 



X 



Here 



0.85 



Sr. 



2 i A{q) B(K,] 



Cc - Cr 

const, and is from (H) with -q - 



(49) 



Figure 7. Critical values of K\ whore outburst amplitudes re- 
sulting from a nova event just vanish, plotted as function of M2 
calculated with Mi = 1.2 Mq. Shown are the lines for orbital 
separation a (full), orbital period P (dotted) and Roche radius 
/Jr (dashed, same as for mass transfer rate — M2) 



1/q. We emphasize that the ratio XfamlZ-^^ does not depend 
on the mass AMcj ejectedper outburst, rather only on Ki 
and the mass ratio q. Fig. H shows Xfaml/A" as function of 
1/q for various values of Ki. Obviously FAML amplification 
is always small for large q (small M2/M1). For a given Ki 
it increases with decreasing q and formally goes to infinity 
where the system approaches thermal instability (^c = Cr)- 
The full line in Fig. ^ indicates the upper limit for 
X-pAMi^l X if the specific angular momentum the ejected en- 
velope can carry is limited by jcj ^ ji + {R2/ a)a^ijj (see 
Sect. 3.1). 



4 NUMERICAL EXAMPLES 

Here we apply the FAML description introduced in the pre- 
vious section and test the analytical considerations of Sect. 2 
and 3 explicitly with numerical sequences of the secular evo- 
lution of CVs. 

The above FAML model contains Ki , AMig and AMoj 
as free parameters. Unfortunately both observations and 
theoretical models describing the outburst itself yield often 
contradicting and inconsistent estimates for these parame- 
ters. Therefore we restrict the following investigation to a 
simple parameter study to obtain a systematic picture of 
the effect of FAML with a given strength on the secular 
evolution. In particular we present calculations with con- 
stant AMig (= AMoj) and Ki, for various values of AMig 
and Ki. 
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Figure 8. FAML amplification factor Xfaml/^i defined in (^), 
as a function of M2/M1 = l/q, for various values of Ki. The 
vertical lines indicate (thermal) instability. The full line is an 
approximate global upper limit, independent of Ki. 



4.1 Computational technique 

To model the binary evolution numerically we describe the 
secondary star by either full stellar models using Mazzitelli's 
stellar evolution code (e.g. Mazzitelli 1989), or by a simpli- 
fied bipolytrope structure using the generalized bipolytrope 
code (Kolb & Ritter 1992). Both codes have been modified 
to include FAML and to resolve individual nova outbursts. 

The bipolytrope code is about 4 orders of magnitude 
faster than the full code and allows one to compute the 
secular evolution for several Gyr with every single nova out- 
burst fully resolved. Several sequences calculated with the 
full code serve as an independent check of the simplified de- 
scription. Despite a careful calibration of the bipolytrope 
model to full stellar models, allowing a quantitative match 
of the results to usually better than 10%, there are well- 
known limitations of the simplified description (see Kolb & 
Ritter 1992), e.g. the increasing deviation of i^ad from the 
values of actual stars for M2 ^ 0.6 Mq. But these effects are 
negligible for the purpose of this paper. 

In all the examples shown below we comp ute magnetic 
braking according to Verbunt & Zwaan (1981) with the cal- 
ibration parameter set to unity. 

We have performed calculations with the following two 
modes: 



Figure 9. Section of a fully resolved evolutionary calculation 



(i) FULLY RESOLVED: Mass transfer is conservative until 
the mass accreted on the WD exceeds the ignition mass. In 
the time step immediately thereafter the ejecta mass and 
the angular momentum ( |39[ ) is removed from the system. 
Thus fully resolved evolutions are discontinuous and show 
repeating nova cycles. 

(ii) CONTINUOUS WIND AVERAGE: Mass loss from the sys- 
tem is continuous at a rate equal to the transfer rate, carry- 



with nova outbursts after AM;™ 



IO-^Mq with 



Ml = 1.2 Mq, M2 ~ 0.6 Mq and weak FAML (Ki = 0.5). 
Crosses indicate time intervals of 1000 years, the filled diamond 
the post-outburst model from which initial values for the ana- 
lytical approximation (dotted line) were taken. The additional 
horizontal lines are explained in the text. 
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Figure 10. As Fig.H, but with strong FAML {Ki = 0.1). 



i ng the increased specific angular momentum according to 
([l4[); hence the secular evolution is contiimous. 
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4.2 Detailed examples for individual outbursts 



We generated short stretches of fully-resolved evolutionary 
sequences with full stellar models by switching into the fully 
resolved mode in the middle of a continuous wind average 
calculation for a typical CV above the period gap. The mode 
switch was made from a model well-established on the uni- 
form evolutionary track characterized by (p^). 

Fig. ^ shows for the case of weak FAML with Ki = 0.5 
a 'sawtooth' (or 'shark-fin') like modulation of the mass 
transfer rate over the orbital period (full line). The system 
evolves from right to left; the crosses mark intervals of 1000 
yr. As expected from Fig. |^ the outburst results in a longer 
period and a lower mass transfer rate with this choice of 
parameters. In the case of strong FAML with K\ — 0.1 
(Fig. |l^) the outburst leads to a shorter period and a higher 
mass transfer rate. 

The approximate time evolution according to Eq. (|l5|), 
with A and B'^ evaluated from the full sequence at the post- 
outburst position, marked with a filled diamond, is shown as 
a dotted curve and matches the full sequence very well. The 
dash-dotted line in Figs. ^ & |l^ indicates the level of the 
'local' stationary conservative mass transfer rate as given by 
Eq. ( p^ ) with — Note that there is a difi'erence be- 
tween this value and the mass transfer rate of a completely 
conservative evolution with the same initial values (but no 
outbursts) as the system here does not evolve conservatively 
on average but loses mass together with angular momentum. 
The averaged stationary mass transfer rate which includes 
FAML with ("r according to Eq. ( [2l| ) is shown as a dashed 
line. This value is identical to the mass transfer rate cal- 
culated in the averaged mode and to the time average of 
the full evolution. The system approaches the conservative 
value between outbursts and oscillates around the average 
transfer rate. 

Note that if, beginning with negligible FAML, the 
strength of FAML is continuously increased, the outburst 
amplitudes will decrease until the stationary (local) conser- 
vative transfer rate and the continuous wind average (which 
grows with FAML) are equal. Further increase of FAML will 
now lead to growing amplitudes, but with opposite signs. 



4.3 Long— term evolution with FAML 

Fig. ^ compares the continuous wind average evolution for 
a reference system (Afi — 1.2 Mq; M2 = 1.0 A4q at turn- 
on; 7; = 0) with strong FAML and without FAML, com- 
puted with full stellar models. As expected from Fig. ^ the 
mass transfer rate with FAML is significantly larger than 
in the case without only at long periods, where the sys- 
tem is close to (but not beyond) thermal instability. Due to 
the higher mass transfer rate the secondary in the FAML 
sequence is driven further out of thermal equlibrium and 
therfore larger at the upper edge of the period gap. At the 
same time the secondary's mass is smaller upon entering 
the detached phase. Both effects, well-known from system- 
atic studies of CV evolution (e.g. Kolb & Ritter 1992), cause 
a wider period gap by both increasing the period at the up- 
per edge and decreasing the period at the lower edge of the 
gap. 

The analytical considerations of Sect. 3.3 are confirmed 
by a parameter study performed with the bipolytrope code. 



u 

o 



o 




Figure 11. Averaged mass transer rate X versus orbital period 
P for a CV with Afi = 1.2 Mq and M2 = 1.0 Mq at turn- 
on without FAML (full line) and with strong FAML (Ki = 0.1, 
corresponding to t)exp — 40 km/sec; dotted line). Computed with 
full stellar models. 



lary 

picted in Fig. Il3 summarize the results. The reference sys- 
tem is shown in the left column, a less massive system in the 
right column. The strength of FAML increases from top to 
bottom, with consequences for the outburst amplitudes as 
discussed in Sect. 3.2 and in the previous subsection. Note 
that the 'sawtooths' here do not represent single outburst 
cycles, but arise because we connect the pre- and post- 
outburst values for every n*'' (n — 1000 or 100) outburst 
only. The advantage of this display method is that both 
magnitude and sign of the outburst amplitude (the verti- 
cal flanks) can be followed through the evolution. In real- 
ity the system undergoes many individual cycles like those 
in Figs. ^ & ^ between subsequent outbursts depicted in 
Fig. Ui. The two evolutionary tracks plotted in each panel 
differ only in the ignition mass AMig (here — AMoj). The 
thick line corresponds to AA/ig = 5 x 10~® Mq, the thin line 
to AAfig = 10~* Mq . In the scale used the low ignition mass 
sequences are practically indistinguishable from the contin- 
uous wind average evolution (remember that both oscillate 
around the same average value, cf. Sect. 4.5). Inspection of 
Fig. 1^ shows that the growth of amplitudes towards short- 
est periods is mainly due to the increase of q. Only in the 



high FAML (low Uoxp) case, i.e. in the very steep region of 
Eq. (^l|), there is an additional contribution from the in- 
creasing orbital velocity (hence decreasing Ki) at short P. 



4.4 FAML— induced instability 

To illustrate and test the analytical mass transfer stability 
considerations of Sects. 2.3 and 2.4 explicitly we consider the 
secular evolution of a system which runs into an instability. 
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Figure 12. Set of fully resolved secular evolutions of a CV witli Mi = 1.2 Mq (left) and Mi = 0.65 Mq (right), computed with the 
bipolytrope code. From top to bottom the strength of FAML increases as labelled. The evolutionary tracks connect the location immedi- 
ately before and after an outburst, hut not for subsequent outbursts! Rather only every every lOOO"" outburst (thick line, corresponding 
to AMig = 5 10~^ Mq), or every 100*'' outburst (thin line, corresponding to AMig = 10~^Mq) is shown. The initial secondary mass at 
turn-on was M2 = 1-0 Mq for the top three sequences on the left, and 0.4 Mq on the right. For the sequences in the lowest panel the 
initial secondary mass was 0.85 Mq and 0.25 Mq, respectively. 



i.e. violates the formal stability criterion (^2|) at some point 
of the evolution. 

To establish this situation we artificially increase the 
strength of FAML along a fully resolved evolutionary se- 
quence (computed with the bipoltrope code, with constant 
AMig = AMej = 1O~*M0), starting from a small value at 
the onset of mass transfer. As a consequence ("j^ becomes 
larger than C^ad at some intermediate secondary mass. In 
particular, we adopt 



km/s 



180 



Ma 
M„ 



50 



(50) 



as the functional form of the envelope expansion velocity. 
(M2 will be large enough to avoid negative Wexp). The cor- 
responding with lypAML from (^ becomes larger than 
("ad at M2 ~ 0.55 Mq, and the mass transfer rate in- 



deed begins to grow without limit at about this mass well 
above the period gap (Fig. |l|). A closer look at the fi- 
nal phase immediately before the runaway (Fig. shows 
that the system evolves beyond the formally unstable point 
apparently unaffected and enters the runaway mass trans- 
fer phase only ~ lO^yr later. In fact this is just as we 
would expect if the evolution were replaced by the corre- 
sponding continuous wind average. In this case (|l^) with 
B — describes the variation of the (average) mass trans- 
fer rate X. The characteristic timescale on which X changes 
is X/X = 1/{A — B^X), which becomes very large when 
B^ ~ 0, i.e. around ("ad = Cr the instability point; note 
that A is also small at this point, see (|3c|). Proceeding fur- 
ther increases the strength of FAML even more and makes 
B^ more negative, hence the timescale progressively 
shorter. The time delay until the runaway finally begins is 
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Figure 13. Fully resolved secular evolution of a CV with M\ = 
0.8 M0 and M2 = 0.7 Mq at turn-on, computed with the bipoly- 
trope code including FAML based on ( |50| ) and AAfig = AAf^j = 
10-'»Mq. For log{P/h) > 0.67 every 100^" outburst is shown, for 
log(P/h) < 0.67 every outburst. 



Figure 14. Evolution of the mass transfer rate X with time dur- 
ing the final phase (approximately 2 million years) before the 
runaway (dotted, filling the grey— shaded area). Overplotted in 
full linestyle is the stationary value Xf = A./B^ for a con- 
tinuous FAML-increased mass loss, indicating an instability at 
~ 12.1 Myr. 



essentially determined by the crossing angle between the 
functions Cad(Af2) and Cr(-^2) and the value of X at this 
point. 

Figure |l^ shows the growth function T defined in ( |27| ) 
for selected models immediately before the runaway for the 
evolutionary sequence depicted in Fig. |l^. The different 
curves are computed with data from the run taken at the 
times as labelled (cf. Fig. As discussed in Sect. 2.4 the 
change of J- from a stable curve with negative slope to an 
unstable curve with positive slope occurs when (^ad = Cr- 

The sequence shown in Figs p^ - p^ is typical of the be- 
haviour of evolutionary sequences with nova cycles and very 
strong FAML. It convincingly demonstrates that the FAML- 
amplified istropic wind average properly describes the effect 
of a sequence of nova outbursts, even in the most extreme 
case when the system approaches a FAML-induced insta- 
bility. A more systematic investigation of the significance 
of this instability for the secular evolution of CVs will be 
published elsewhere (Kolb et al., in preparation). 



4.5 Evolution with FAML parameters from 
theoretical nova models 

The FAML parameters K\, AAfoj and AAfig are certainly 
not constant along the evolution but depend on the actual 
state of the outbursting system. The three governing pa- 
rameters which crucially determine the outburst character- 
istics of thermonuclear runaway (TNR) models for classical 
novae are the WD mass M\, the (mean) accretion rate X 
and the WD temperature Ti (e.g. Shara 1989). Both the 
mass and the temperature of the WD change only slowly 
during the binary evolution, thus the dominant dependence 
of Ucxp is from the (mean) mass transfer rate. As FAML 



itself drives the mass transfer this could lead to an interest- 
ing feedback and possibly self-amplification in the system, 
a question considered in a separate paper (Kolb et al., in 
preparation). 

Although published TNR models mostly give values for 
AAfig and estimates for AA/oj, few or no data are available 
on Ki. Unfortunately this is also true for the most com- 
plete and consistent set of published nova models (Prialnik 
& Kovetz 1995) which we chose to use as input to the FAML 
description derived in Sect. |j.l| . 

These authors tabulate only the average expansion ve- 
locity Uav , the time average of the maximum velocity in the 
flow during the whole mass loss phase. In the absence of 
more relevant velocity data we simply set iiexp = Vav, but 
we are aware that this choice represents an upper limit for 
«cxp (which has to be taken at the secondary's position). 
Two exemplary relations Vcyip{X) are shown in Fig. ^ for a 
low and a high mass WD, complemented by large values (10^ 
km/sec, hence negligible FAML) when no envelope ejection 
was found in the hydrodynamical simulations. 

The evolutionary sequences obtained with M\ = 1.2Mq 
(using the velocity data for nova models with a 1.25 Afp 
WD) and with Mi = 0.65 A/q are shown in Figs. |l^ & |l^ 
The changes of Afi during the whole evolution are very small 
(77 ~ 0). As in Fig. ^the 'sawtooths' result from connecting 
only every n"^ outburst. For comparison a standard contin- 
uous wind average evolution without FAML but otherwise 
identical parameters has been overlaid in thick linestyle. As 
expected from the high expansion velocities FAML has very 
little effect on the global evolution. By choice we underesti- 
mated the strength of FAML in the sequences shown above. 
But although the expansion velocity at the location of the 
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Figure 15. Growth function T (cf. ( |28[ ) and Fig. |^) calculated 
with data talcen from the evolution shown in the previous figures 
for different times as labelled. 



secondary is ind eed s maller than Wav in the models by Pri- 
alnik & Kovetz (1995) they currently appear to be too large 
to cause a FAML effect with significance for the long-term 
evolution of CVs (cf. also discussion in next section). 

Perhaps the most interesting feature of the sequences 
presented in this paragraph is the large outburst amplitude 
of the mass transfer rate in low-mass WD (O.GSM©) systems 
below the period gap. The reason for this lies mainly in the 



fact that AA/cj and AAfig to a first approximation scale like 
Rf/Mi (Fujimoto 1982, Ri is the WD radius), i.e. increase 
with decreasing WD mass (and to a lesser extent, increase 
with decreasing X). Moreover the slow but continuous de- 
crease of the WD mass leads to a further increase of the 
outburst amplitudes towards shorter orbital period in this 
sequence (Fig. ^sj). Note that this effect adds to the ones 
observed before (lowermost panels of Fig. ^2|) which were 
caused by the change of q due to decreasing M2, and to a 
lesser extent to the slowly increasing orbital velocity Vscc, 
mimicing increasing FAML strength. 

Such amplitudes account for a non-negligible time in- 
terval with a secular mean mass transfer rate significantly 
below the continuous wind average transfer rate, i.e. below 
~ 5 X lO~^^M0/yr, the typical value for non-degenerate 
CVs below the period gap driven by gravitational radiation 
alone. As an example Fig. ^ plots for a model at P = 1.23 h 
from the sequence shown in Fig. |l^ the time spent with a 
mass transfer rate below a value X, as a function of X. One 
nova cycle lasts ~ 6.5 Myr, and the mass transfer rate re- 
mains for 0.5 Myr after the outburst below lO~^^Af0/yr. 
Intrinsically, these CVs with low-mass and moderate-mass 
WDs form the vast majority (e.g. de Kool 1992, Politano 
1996). Thus, ignoring selection effects, we would expect to 
observe 1 out of 10-20 short-period CVs with small sub-GR 
driven mass transfer rate. 

This seems to offer an alternative explanation for the 



Figure 16. Averaged maximal expansion velo citv o f the ejecta 
taken from nova models of Prialnik & Kovetz (1995), with Ti = 
10^ K and Mi = 0.65 Mq (full line) or Mi = 1.25 Mq (dot- 
ted line). Diamonds mark real data points. A maximum value of 
10^ km/sec has been inserted if no ejection occurred to mimic 
negUgible FAML. 
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Figure 17. Secular evolution of a CV with fully resolved nova 
cycles and Mi = 1.2 Mq, M2 = 1.0 Mq at turn-on (thin hne), 
computed with the bipolytrope code. Ignition mass , ejec tion mass 
and velocities are taken from Prialnik & Kovetz (1995). Similar 
to Fig. ^ only every lOO*'' outburst is shown, and straight lines 
connect post— and pre— outburst points. For comparison the con- 
tinuous wind average evolution without FAML (thick line) is also 
shown. At the period minimum the WD mass has decreased to 
Ml ^ 1.1 Mq. 
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Figure 18. As Fig. |l7[ but with Mi = 0.65 Mq, Nh.i = 0.4 Mq, 
and every 20"^ outburst shown. At the period minimum the WD 
mass has decreased to Mi ~ 0.635 Mq. 
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Figure 19. A typical inter-outburst low— mass WD CV below 
the period gap (Mi = 0.6382 Mq, M2 ~ 0.1 Mq, AM^j = 
AMig 
alnik 



2.9 10 MfTi) according to nova parameters from Pri- 



Kovetz (1995). Shown is the time spent with a mass 



transfer rate below a value X, as a function of X (in Mq/jt). Di- 
amonds indicate the immediate post- or pre— outburst state. The 
dashed line indicates the continuous wind average mass transfer 
rate. 



low mass transfer rate Sproats et al. ( 199(^ ) claim to find 
observationally in so-called tremendous out burst amplitude 
dwarf novae (TOADs) which Howell et al. (1997) interpret 
as post minimum period systems. As the secular mean mass 
transfer rate is predicted to drop substantially (j5 10"^'^ — 
lO~^^M0/yr) CVs which have evolved past the minimum 
period would be even fainter than the above post-outburst 
CVs, and it is unclear if they are detectable at all. 

We n ote that at least 4 objects in the hst of Sproats et 
al. (1996) have been labelled as novae (cf. Duerbeck 1987): 



AL Com (orbital period 0.0566 d, outburst suspected 1961), 
VY Aqr (0.0635 d, 1907), RZ Leo (0.0708 d, 1918) and SS 
LMi (no period known; 1980). With exception of SS LMi 
they are now firmly considered to be dwarf novae. This does 
not exclude a priori that the first recorded outburst (e.g. 
in 1907 for VY Aqr) or even earlier ones were actually no- 
vae. Hence it is at least conceivable that the present low 
mass transfer rate and the TOAD characteristics might be 
a consequence of this last nova event. The three systems 
with known orbital period are lo cated at the extreme end 
of Fig. 4 of Sproats et al. (1996), showing large outburst 
amplitude and low quiescence magnitude. Moreover, among 
the open square symbols indicating DNe in their Fig. 3 they 
are also those with lowest period and absolute magnitude 
in quiescence (mass transfer rate -ti 10~'^^MQ/yr). The fact 
that the claimed deviation from the secular mean magnitude 
is shrinking with the time elapsed since the potential nova 
outburst (1961 - 1918 - 1907) is probably just a coincidence. 
We caution that it is difficult to decide if the set of 



models by Prialnik & Kovetz (1995) properly describes clas- 
sical nova outbursts on WDs in CVs below the period gap. 
Observationally, only 4 out of 28 classical novae with deter- 
mined orbital period are below the gap (e.g. Ritter & Kolb 
1995). This is in confiict with standard population models 
of CVs if the Prialnik & Kovetz ( |l995| ) ignition masses are 
used to predict an observable period distribution for novae 
(Kolb 1995). 



5 DISCUSSION AND CONCLUSIONS 

In this paper we considered the effects of nova outbursts on 
the secular evolution of CVs. As a result of these outbursts 
the secular mean mass transfer rate and the orbital period 
are not continuous functions of time but change essentially 
discontinuously with every nova outburst by an amount pro- 
portional to the ejected envelope mass. In addition, energy 
and angular momentum can be removed from the orbit due 
to dynamical friction of the secondary orbiting in the ex- 
panding nova envelope. 

The discontinuous evolution with a given strength of 
frictional angular momentum loss (FAML) is usually re- 
placed by the corresponding continuous wind average evolu- 
tion, where the mass and angular momentum loss associated 
with a nova outburst is assumed to be distributed over the 
inter-outburst time and to form an isotropic wind from the 
white dwarf. We showed analytically that the well-known 
mass transfer stability criterion for the latter case can also 
be derived from a proper analysis of the real, discontinuous 
process, for an arbitrary strength of FAML-amplification. 

We specified a quantitative model for FAML within the 
framework of Bondi-Hoyle accretion following Livio et al. 
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(1991). In this model the strength of FAML depends cru- 
cially on the expansion velocity Uexp of the envelope at the 
location of the secondary, being the stronger the smaller iiexp 
is. We expect that the resulting simple one-parameter de- 
scription properly describes the order of magnitude of the 
FAML effect and, more importantly, the differential depen- 
dences on fundamental binary parameters. Hence although 
it is a useful way to study the potential influence of FAML 
systematically, it certainly cannot replace a detailed mod- 
elling of the frictional processes. 

Calculations of the long-term evolution of CVs verified 
the validity of the replacement of the discontinuous sequence 
of nova cycles with the continuous wind average, even for 
situations close to mass transfer instability, whatever the 
strength of FAML. The mass transfer rate in the continu- 
ous wind average evolution is FAML-amplified, i.e. by the 
factor ( |4E| ) larger than the transfer rate driven by systemic 
angular momentum losses alone. For a given CV this factor 
is determined by K\ alone, i.e. independent of the ejection 
mass. We emphasize that FAML only amplifies the trans- 
fer rate caused by systemic losses, it does not add to them. 
Hence FAML is a particular example of consequential an- 
gular moment um lo sses (CAML) investigated in detail by 
King & Kolb ( |1995| ). 

In general, the FAML amplification factor turns out 
to be large only when the envelope expansion is very slow 
{Ki ^ 0.5, i.e. iiexp ^ 200 km/s) and when the system is al- 
ready close to thermal mass transfer instability (Fig.^ . This 
latter condition means that e.g. for an evolution with strong 
FAML Ki = 0.1 = const, the averaged mass transfer rate is 
significantly affected only at long orbital period, P ^ 5 h. 

The magnitude and direction of the outburst ampli- 
tudes of the mass transfer rate X and the orbital period 
P depend on both the ejection mass AAfcj and the FAML 
parameter Ki . For weak (or negligible) FAML the outbursts 
are towards lower mass transfer rate and longer orbital pe- 
riod, for strong FAML towards larger X and shorter P. 
There is an intermediate regime where the outburst am- 
plitudes essentially disappear. 

Theoretical models for nova outbursts generally find 
larger expansion velocities (e.g. Prialnik 1986, Prialnik & 
Kovetz 1995), typically Ki J> 1. This is certainly true for 
the terminal velocities, but in more recent models probably 
also for the crucial velocity at the s econd ary's location, i.e. 
closer to the WD. Kato & Hachisu (1994) argue that in the 



wind mass loss phase of a nova outburst the main accelera- 
tion (at about the sonic point) takes place at a temperature 
where the opacity has a maximum. Thus the introduction 
of the OPAL opacities had the effect of moving this point 
closer to the WD. A comparison of th e radial velocity pro- 
files found by Kato & Hachisu ( 1994 ) with those obtained 
previously (e.g. Prialnik 1986, Kato 1983) confirms this. As 
a result, in the newer models the velocities at the location 
of the secondary are already quite close to their terminal 
values. This explains why Kato & Hachisu (1994) find only 
marginal effects from dynamical friction. 

This seems to suggest that the overall infiuence of 
FAML on the long-term evolution of CVs is small. However, 
in view of the considerable simplifications of our FAML de- 
scription and the uncertainties of theoretical TNR models 
for nova outbursts, it is worthwhile to investigate sytemati- 
cally mass transfer stability with FAML, and to consider the 



role of a feedback between outburst characteristics (hence 
FAML strength) and the mean mass transfer rate prior to 
the outburst. We will study this in a forthcoming paper 
(Kolb et al., in preparation). 

To illustrate further the effect FAML might have on the 
secular evolution of CVs we have calculated evolutionary 
sequences with envelope expansion velocities, ignition and 
ejection masses taken fr om t he extended set of nova models 
by Prialnik & Kovetz (1995). As expected, the continuous 



wind average evolution hardly differs from the standard CV 
evolution without FAML. As a consequence of the large ejec- 
tion mass (several ~ lO~*M0/yr) the outburst amplitudes 
become very large (a factor ^ 10) below the period gap for 
intermediate mass WDs (Mi ~ O.6M0) — even more so as 
the outburst-induced decrease of the mass transfer is largest 
if FAML vanishes. 

Such systems have a mass transfer rate less than 
lO~^^M0/yr for ~ 0.5 Myr after the outburst, and could 
account for intrinsically faint CVs below the period gap. We 
speculated (Sect. 4.5) if TOADs (e.g. Sproats et al. 1996) 
could represent such systems. 

We finally note that the apparent scatter in observation- 
ally derived values for the mass transfer rate of CVs with 
comparabl e orbi tal period, well-known since the review of 
Patterson (1984), is unlikely to be due to outburst ampli- 
tudes. First, we expect from Fig. ^ that these amplitudes 
are small or negligible above the period gap, and second 
the systems spend most of the time close to the continu- 
ous wind average mass transfer rate (cf. Fig. [l9| ). A more 
promising explanation for this scatter assumes mass trans- 
fer cycles which could be irradiation-induced (e.g. King et 
aL 1995, 1996). 
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(/32 + 2)(7-l) 



which with (0, (|l|) and AM2 = aAAfoj - (7 - l)AMig 
gives the generaUzed 
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Six parameters describe the system: AAfig and AA/oj 
(ideally taken from nova models), ^taml and a (requiring 
the specification of a particular FAML model), and 7 (or 
vT), (characterizing the wind loss during accretion). 

The ratio of mass retained by the WD to that lost by 
the secondary 

/AMi\ _ AAfcj-AAfig 



VAA/a 



QAA^cj - (7 - l)AA^ig ' 



(A4) 



which relates AAfoj and AA/ig similar to the simple case rj" 
in the main body of the paper, now depends on a and 7. 
Using a weighted u to describe the average specific angular 
momentum gives 



(a - l)AAfej (i + i^faml) - (7 - 2)AAfig^'' 
{a - l)AA^oj - (7 - 2)AA/ig 



(A5) 



and inserting these values into Eq. (g) would also have di- 
rectly led to Eq. (O). 



APPENDIX A: GENERALIZATION OF 



Starting from the fundamental Eq. (j^ one can derive a more 
general expression for the relative change of _Rr during both 
the outburst and inter-outburst phase. Specifically we allow 
a small fraction a of the mass ejected during the outburst to 
accrete on to the secondary, i.e. AA/2 ~ ol AAfcj . This gives 



/ AitR \ 
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-; (1 H h 2 ^FAMLJ 

1 + <? q 

instead of (p7|), which is obtained in the case a — 0. 

Additionally we allow a wind (or other) mass loss during 
the inter-outburst accretion phase. Thus we use the general 
^R, from (^ rather than the conservative one to general- 
ize (^^. This results in two additional parameters, rj™ and 
, describing the system's "wind" mass and corresponding 
angular momentum loss during the inter-outburst phase. 
Hence we get (using 7 = {rj"' + l)/??"') 
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